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Abstract 

We review the connection between inflationary models and obser- 
vations and concentrate to describe models based on softly broken 
supersymmetry, in particular running mass models, and their predic- 
tions. We then present a fit of the spectral index of the curvature 
perturbation, assuming a flat ACDM cosmology. 

1 Introduction 

An epoch of inflationary expansion of the Universe is necessary to solve some 
of the problems of the standard Big Bang cosmology, like the large scale ho- 
mogeneity and isotropy, flatness and unwanted relics problems. Moreover it 
has been also found that slow roll inflation can generate the small scale struc- 
ture in the Universe through a primordial gaussian curvature perturbation, 
originated from the quantum fluctuations of the inflaton field. 

In general a model of inflation consists in a scalar potential V(4>) for the 
inflaton field <fi satisfying slow roll conditions^] 
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V = M p[y) «! ( 2 ) 

where the prime denotes derivative with respect to the field <f>. 

The point of contact between observation and models of inflation is the 
spectrum Vn(k) of the curvature perturbation, which, in the slow roll ap- 
proximation 3H(p ~ —V, is given in terms of the inflaton potential V{<p) 

1 As usual, in the following Mp = 2.4 x 10 18 GeV is the Planck mass, a is the scale 
factor and H = a/a is the Hubble parameter, and k/a is the wavenumber. 
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by 

where the potential and its derivatives are evaluated at the epoch of horizon 
exit k = aH. To work out the value of 4> at this epoch one uses the relation 

Hk cnd /k) = N(k) = Mp 2 j* (V/V)d4, (4) 

where N(k) is actually the number of e-folds from horizon exit of the scale 
k to the end of slow-roll inflation. At the scale explored by the COBE mea- 
surement of the cosmic microwave background (cmb) anisotropy, N{kcoBE) 
depends on the expansion of the Universe after inflation in the manner spec- 
ified by: 

Ncobe * 60 - ln(10 16 GeV/U 1 / 4 ) - ± HV^/T^) . (5) 

In this expression, T Te ^ is the reheat temperature, and instant reheating is 
assumed. 

1/2 

Given the above relations, the observed large-scale normalization — 
1CP 5 provides a strong constraint on models of inflation. But here we are in- 
terested in the scale-dependence of the spectrum, defined by the, in general, 
scale-dependent spectral index n; 

dlnVTi(k) . . 

According to most inflationary models, n has negligible variation on cosmo- 
logical scales so that Vn(k) oc fc n_1 , but we shall also discuss an interesting 
class of models giving a different scale-dependence. 
From (|) and @, 

n-1 = 2M 2 P (V"/V) - 3Mj,{V' /V) 2 , (7) 

and in almost all models of inflation, (0) is well approximated by 

n - 1 = 2M 2 P {V"/V) . (8) 

We see that the spectral index measures the shape of the inflaton potential 
V (4>), being independent of its overall normalization. For this reason, it is 
a powerful discriminator between models of inflation. 
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2 Supersymmetric models of hybrid inflation 



A very successful and interesting type of inflationary model is hybrid infla- 
tion, proposed by Linde ||. In this case the potential depends on two fields, 
the inflaton and another field, which triggers a phase transition and the end 
of slow roll inflation and acquires a non vanishing v.e.v. at the end. The 
basic potential is of the form: 

v = vM + \{<P 2 - m i)^ + ---- ( 9 ) 

so, for <j) 3> m^p, the hybrid field tp has vanishing v.e.v. and inflation can take 
place driven by V^, but when the inflaton becomes smaller, -ip is displaced 
from the origin and its v.e.v. stops slow roll inflation. Then a phase of 
oscillation of the two fields around the true minimum of the potential occurs 
and the preheating/reheating process. 

We see then that the characteristic of hybrid inflation is that two differ- 
ent part of the potential are responsible for the inflationary phase and for 
the following stage of reheating. So the power spectrum of the curvature 
perturbations depends only on the inflationary potential and the critical 
value (p c = m^p, which corresponds to the end of slow roll inflation. 

Supersymmetric versions of hybrid inflation have been studied by many 
authors (see [pj and references therein) and are based on the superpotential 

W = A$(Af 2 - f^), (10) 

where capital Greek letter denote the superfields corresponding to the lower 
case Greek letter scalar field. 

In the limit of global super symmetry, such potential is perfectly flat 
with respect to the inflaton field and the small slope needed for slow roll is 
generated by supersymmetry breaking effects. Depending on the different 
mechanism considered, different scenarios are possible, as described in Table 
H|, together with the prediction for the spectral index and its derivative. 

Note that the COBE normalization constraints the overall scale Vo, while 
the spectral index gives direct informations on the supersymmetry breaking 
effect which dominates the inflaton potential. 
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Table 1: Various models of hybrid inflation and their prediction for the spec- 
tral index as a function of N = —log(k/k en d) and its derivative, assuming 
small field values. In the case of the simple linear term, n — 1 is given by the 
full expression (0) since rj vanishes, but note that in supergravity generally, 
other contribution to the spectral index coming from higher order terms are 
usually present and can be larger than the one listed here Q. 

Vj>((f))/Vo n—1 dn/d\og{k) Origin of the 

slope 



& log — ^ , -\ 2V 1 W for 

N+ 2 Z 4 'A ) spont. broken 



4^ &V Q J N+' 27r2 ' p 



susy 

1±\^4> 2 Susy breaking 

mass 



l±Pjfe -3/? 2 Susy breaking 

linear term 

1 ± t*t .J 2 . J 12 r i Sugra quartic 



^ N [ ^ T N ) term 
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3 Running mass models 



3.1 The potential 

Let us consider now a specific type of models of hybrid inflation, the models 
with a running inflaton massj|, [|, ^, In these models, based on softly 
broken super symmetry^, one-loop corrections to the tree-level potential are 
taken into account by evaluating the inflaton mass-squared m 2 (ln(Q)) at 
the renormalization scale Q ~ eft, as long as <j> greater than all other relevant 
scales. 

V = V Q + ^m 2 (\n(Q))cf ) 2 + --- . (11) 

The dependence on Q is only logarithmic due to super symmetry, which 
ensures the cancelation of the quadratic divergences; anyway since super- 
symmetry is broken by the soft terms, a logarithmic divergence survives. 

Over any small range of 4>, it is a good approximation to take the running 
mass to be a linear function of In <f> and then one reproduces explicitly the 
expression of the tree potential plus loop correction: 



V = V + im 2 (ln(Q))0 2 - I c (fn(Q))^ 2 ]n(</>/Q) + ■■■ . ( 12 ) 



It has been shown ||] that the linear approximation is very good over 
the range of <f) corresponding to horizon exit for scales between kcoBE and 
8fo _1 Mpc. We shall want to estimate the reionization epoch, which involves a 
scale of order k~^ ion ~ 10~ 2 Mpc (enclosing the relevant mass of order 10 6 ©). 
Since only a crude estimate of the reionization epoch is needed, we shall 
assume that the linear approximation is adequate down to this 'reionization 
scale'. 



A potential of the type (12) gives rise to four different models of inflation, 
depending on the sign of c and the direction of rolling of the inflaton field, 
towards or away from the origin^. 

2 We denote as spontaneously broken supersymmetry any theory where supersymmetry 
breaking preserves the super-trace relation StrM = 0, while in softly broken supersym- 
metry one has Str.A/! 2 7^ 0; e.g. the displacement of the inflaton field from the minimum 
of the potential causes a spontaneous breaking in the inflaton sector (which can give a 
soft breaking in the visible sector), while spontaneous or explicit supersymmetry breaking 
in another (hidden) sector is needed to have soft breaking in the inflaton sector. Note 
that in the first case, supersymmetry is restored once the inflaton field settles in the true 
minimum, while in the second case the breaking persists all the time, even if it can be 
affected by the dynamics during inflation. 

3 We follow the labeling introduced in M. 
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The value of c is given by the well-known Renormalization Group Equa- 
tion for the inflaton mass, and it depends on the gauge and Yukawa cou- 
plings, a and A respectively, of the inflaton field: 

V c dm 2 2C ^ 2 D 2 2 

~Mj = dT = ~ am ~ I6^ |A| m '«* • (13) 

Here, C is a positive group-theoretic number of order 1, fh is the supersym- 
metry breaking gaugino mass, D is a positive constant counting the number 
of scalar particles interacting with the inflaton through A and mf oop is their 
common susy breaking mass-squared. Negligible supersymmetry breaking 
trilinear coupling is assumed in the above expression. 

We see that the running is then directly connected to the couplings of 
the inflaton field and the supersymmetry breaking mass spectrum of our 
model and so to complete our estimate of c, we need the gaugino or scalar 
mass. A very minimal and traditional hypothesis is that soft supersymmetry 
breaking is gravity-mediated and that the scale of susy breaking during 
inflation V^ A coincides with the scale of supersymmetry breaking in the 



present vacuum, Ms = \F, where F is the auxiliary field responsible for 
spontaneous supersymmetry breaking in the hidden sector. 

With gravity-mediated susy breaking, typical values of the masses are 
fh 2 ~ \m 2 oop \ ~ Vo/Mp, which makes |c| of order of the coupling strength a 
or |A| 2 . 

At least in the case of dominance of the gauge coupling, one then expects 
a small, but non-negligible value for c 

|c| ~ 10 _1 to 10~ 2 . (14) 

In special versions of gravity-mediated susy breaking, the masses could be 
much smaller, leading to |c| <C 1. In that case, the mass would hardly run, 
and the spectral index would be practically scale- independent. With gauge- 
mediated susy breaking, the masses could be much bigger; this would not 
lead to a model of inflation (unless the coupling is very suppressed) because 
it would not satisfy the slow-roll requirement |c| < 1. 

So we see that both from the theoretical and from the observational 
perspective, the interesting parameter region for the running mass models 
is given by ([ill). 

3.2 The spectrum and the spectral index 

Using (|j) we find 

se cAN(k) = cln (^/0) ( 15 ) 
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AN(k) = NcoBE-N(k) = ln(k/k C oBE), (16) 



where s is an integration constant, which absorbs the dependence on the 
end of inflation. The spectral index is then given, using (||), by 

T ^l_l =se ^N {k )_ c _ 

We see then that a relatively strong scale dependence arises for n(k), de- 
pending on the magnitude of the constants s and c. To satisfy the slow-roll 
conditions (1), both c and s must be smaller than 1 in magnitude. 

In order to evaluate the power spectrum later on (see ( |27l) in the following 
section) , we also need the variation of 5h which comes from integrating this 
expression, 



S H (k) 

— exp 



£ cAN 



, i • ■ / (18) 

0H{kcOBE) 

We are mostly interested in cosmological scales between kcoBE and kg, 
corresponding to < AN < 4. In this range the scale-dependence of n is 
approximately linear (taking |c| < 1) and the variation An = n% — ncoBE 
is given approximately by 

din 

An~4— — -~8sc. (19) 
dmk 

4 Fit to a ACDM model 

We will present here a global fit to the cosmological parameters according 



to the procedure described in full detail in |H|. Following [10[ we will take 



into account also the recent measurements of the cmb anisotropy performed 



by the Boomerang [11] and Maxima-1 [12| balloon experiments. 

The observational constraints on the cosmological parameters and on 
the spectral index have been studied by many authors, but in our analysis 
we will use a different treatment of reionization and we will compare the 
present data with the two parameters prediction for the scale dependence of 
the spectral index in running mass models. 

Observations of various types indicate that we live in a low density Uni- 



verse, which is at least approximately flat (see e.g. [11, |12|)- In the interest of 
simplicity we therefore adopt the ACDM model, defined by the requirements 
that the Universe is exactly flat, and that the non-baryonic dark matter is 
cold with negligible interaction. Essentially exact flatness is predicted by 
inflation, unless one invokes a special kind of model, or special initial con- 
ditions. We shall constrain the parameters of the ACDM model, including 
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the spectral index, by performing a least-squares fit to key observational 
quantities. 



4.1 The parameters 

The ACDM model is defined by the spectrum Vn(k) of the primordial cur- 
vature perturbation, and the four parameters that are needed to translate 
this spectrum into spectra for the matter density perturbation and the cmb 
anisotropy. The four parameters are the Hubble constant h (in units of 
100kms _1 Mpc _1 ), the total matter density parameter Qq, the baryon den- 
sity parameter fij, and the reionization redshift zr. As we shall describe, zr 
is estimated by assuming that reionization occurs when some fixed fraction 
/ of the matter collapses. Within the reasonable range / ~ 10 -4 to 1, the 
main results are insensitive to the precise value of /. 

The spectrum is conveniently specified by its value at a scale explored 
by COBE, and the spectral index n(k). We shall consider the usual case 
of a constant spectral index, and the case of running mass models where 
n(k) is given by the two-parameter expression in (|17|). Since VnikcOBE) is 
determined very accurately by the COBE data we fix its value. Excluding 
zr and Vn{kcoBE)i the ACDM model is specified by four parameters in 
the case of a constant spectral index, or by five parameters in the case of 
running mass inflation models. 



4.2 The data 

As described in detail in P, 10 1, we consider a sample of data constituted 



by seven observational quantities. Of these quantities, three are the cosmo- 
logical quantities h, S7o , fig, which we are also taking as free parameters. 
We will assume that, at least at some crude level, we can pretend that the 
errors are all random and uncor related, and perform a least squares fit. 

The other data are taken so to sample different observable scales: one 
one side we use large scale structure observation, on the other measurements 
of the cmb anisotropy. 

So, first of all we consider the rms density perturbation at 8/i _1 Mpc, og, 
measured through the abundance of rich galaxy clusters at redshift z = to 



a few [13 1 and the shape parameter T that specifies the slope of the galaxy 
correlation function on scales of order l/i -1 to 100h~ 1 Mpc [|14|, [l5| 

Secondly we take also two data from the latest measurements of the cmb 
anisotropy, i.e. the height of the first peak at £ ~ 210 — 30 and the ratio 



S 



between the heights of the second and first peak. We consider the average of 



the two experiment [11, 12] and include the calibration error in quadrature. 



We fix the value of the COBE normalization, as described in detail in 

The adopted values and errors are given in Table 2 and 3, with the results 
of the fit for a constant and scale-dependent n. 

4.3 Reionization 

The effect of reionization on the cmb anisotropy is determined by the optical 
depth r. We assume sudden, complete reionization at redshift zr, so that 
the optical depth r is given by |[(| 

r = 0.035^/i ^Jn (l + z R ) 3 + l-n - l) • (20) 

In previous investigations, zr has been regarded as a free parameter, 
usually fixed at zero or some other value. In this investigation, we instead 
estimate Zr, in terms of the parameters that we are varying plus assumed 
astrophysics. Indeed, it is usually supposed that reionization occurs at an 
early epoch, when some fraction / of the matter has collapsed into objects 
with mass very roug hly M = 10q. Estimates of / are in the range |[7|] 

HT 4 - 4 </<l. (21) 
In the case / < 1, the Press-Schechter approximation gives the estimate 

l + ^-^^erfc-H/) (/«!). (22) 

Here a{M) is the present, linearly evolved, rms density contrast with top- 
hat smoothing, and 5 C = 1.7 is the overdensity required for gravitational 
collapse (g is the suppression factor of the linearly evolved density contrast 
at the present epoch, which does not apply at the epoch of reionization. See 
H for details). In the case / ~ 1, one can justify only the rough estimate 

l + z R ~^\ (23) 

not very different from the one that would be obtained by using / = 1 in 
©• 

In our fits, we fix / at different values in the above range, and find 
that the most important results are not very sensitive to / even though the 
corresponding values of zr can be quite different. 
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4.4 The predicted peak height 



The CMBfast package [18] gives Ci, for given values of the parameters with n 
taken to be scale-independent. Following [19], we parameterize the CMBfast 
output at the first peak in the form 



Cpeak — Cpeak [ TrT ) ' ( 24 ) 



/220 V/ 2 



'peak y\Q 

e value of C peak 

where v = 0, and 



where C^J ak is the value of C pea k at the reference value of the parameters, 



v = a n {n-l)+a h ln(/i/0.65) +a ln(O /0.35) + a b /i 2 (O b -0.019)-0.65/(r)T . 

(25) 

The coefficients are a n = 0.88, = —0.37, ao = —0.16, a& = 5.4, 
and C^J ak = 77.5uK. The formula reproduces the CMBfast results within 
10% for a 1-cr variation of the cosmological parameters, h, Qq and f^, and 
ncoBE = 1-0 ± 0.05. 

With the function /(r) set equal to 1, the term — 0.65r is equivalent to 
multiplying C pea k by the usual factor exp(— r). We use the following formula, 
which was obtained by fitting the output of CMBfast, and is accurate to a 
few percent over the interesting range of r; 

/ = l-0.165r/(0.4 + r). (26) 

For the running-mass model, we start with the above estimate for n = 1, 
and rescale it according to the scale dependence of (|^), i.e. 

Cpeak _ 6 H (k(£,n )) ^ 



^("=1) S H (kcOBE(^o)) ' 
peak 

In the case of constant n, this prescription corresponds to the previous one 
with a n = 0.91, in good agreement with the output of CMBfast. 



5 Results 

5.1 Constant n 

For the case of a constant spectral index our results are summarized in Table 
@, together with the data points used. 

The relative high value of the x 2 is mainly due to the second peak con- 
tribution, but note that the discrepancy between the measured and fitted 
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Table 2: Fit of the ACDM model to presently available data, assuming reion- 
ization when a fraction / = 10 -2 ' 2 of matter has collapsed (corresponding 
redshift at best fit is zr = 18). The first two rows show the data points, 
while the result of the least-squares fit for the parameters is given in the 
lines three to five. All uncertainties are at the nominal 1-a level. The total 



2 is 6.3 with three dej 


frees of freedom 










n 


n b h 2 


^0 


h 


r 






r ~i - " d 

°e 


data 


0.019 


0.35 


0.65 


0.23 


0.56 


74.0 \iK 


0.38 


err. 


0.002 


0.075 


0.075 


0.035 


0.059 


5.0 [l K 


0.06 


fit 0.99 


0.021 


0.38 


0.62 


0.19 


0.56 


70.8 \lK 


0.49 


err. 0.05 


0.002 


0.06 


0.05 










x 2 - 


0.9 


0.2 


0.2 


1.3 


0.002 


0.4 


3.3 



values of the peak ratios is within the 2a limit and that the best fit value of 
tl^h? is also within la from the Big Bang Nucleosynthesis measurement. 

We show also in Figure 1, the comparison of the dependence of the 
best fit value of the spectral index on the the reionization redshift with the 
dependence on the reionization fraction /. Note that a strong correlation 
between n and zr is present, while the correlation with / is much weaker, 
and in particular the lower bound on n is practically independent of /. Low 
values of zr and therefore high / have lower x 2 - 

5.2 Running mass models 

We performed the same fit using also the scale dependent spectral index 
given in (|l7|). In this case we obtain an allowed region in the plane s vs c, 
shown in Figure 2. In the same figure it is also displayed the most reasonable 
region of values for the parameters, from theoretical arguments (see ]|] for 
a full discussion). 

In Table 3 are again summarized the results of the fit in this case: we 
have fixed the value of c, so that the number of degrees of freedom is the 
same as before. The value of n in the table, refers to a scale corresponding 
to the COBE measurement. Since c is positive, smaller scales present larger 
spectral index n > ucobe, e.g. n(fcg) = 1.009. Note that the value of \ 2 is 
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Table 3: Fit of the ACDM model to presently available data for the running 
mass models. We have fixed c = 0.1 and so the free parameters are ucobe = 
1 + 2(s — c), and the next three quantities in the Table. Reionization is 
taken to occur when a fraction / = 10~ 2 ' 2 of matter has collapsed, as done 
previously (the corresponding redshift at best fit is zr = 21.) Every quantity 
except ticobe is a data point, with the value and uncertainty listed in the 
first two rows. The result of the least-squares fit is given in the lines three 
to five. All uncertainties are at the nominal 1-cr level. The total \ 2 is 8.4 
with three degrees of freedom 





^■COBE 


n b h 2 
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c\ at 


data 




0.019 


0.35 


0.65 


0.23 


0.56 


74.0 \iK 


0.38 


err. 




0.002 


0.075 


0.075 


0.035 


0.059 


5.0 |lK 


.06 


fit 


0.94 


0.021 


0.40 


0.59 


0.19 


0.53 


67.6 |lK 


0.49 


err. 


0.04 


0.002 


0.05 


0.05 










x 2 




0.9 


0.4 


0.6 


1.2 


0.2 


1.6 


3.5 



in this case larger: the minimum lies on the c = line and favours no scale- 
dependence. Note anyway that along the quasi degenerate s — c direction, 
corresponding to small ucobe, values of |c| < 0.1 — 0.2 are within the 70% 
CL contour. 
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6 Conclusion 



We have presented a fit of the ACDM model to a global data set, assum- 
ing that a gaussian primordial curvature perturbation is the only one. We 
focused on the spectral index n, specifying the shape of the curvature per- 
turbation, considering separately the case of a practically scale-independent 
spectral index, and the scale-dependent spectral index predicted by running 
mass inflation models. In contrast with other groups, we calculate the reion- 
ization epoch within the model on the assumption that it corresponds to the 
epoch when some fraction / of the matter collapses, the results being only 
mildly dependent on / in the reasonable range / > 10~ 4 . 

For the scale-independent case, we obtain a spectral index n = 0.99±0.05 
for an intermediate value of /. The result for all / is given in Figure 1. We 
stress that the best fit result strongly depends on the procedure used to 
treat the reionization redshift zr. 

In the case of running mass models, the scale-dependent spectral index 
depends on parameters s and c, the latter being related to the inflaton 
coupling which produces the running. We have delineated the allowed region 
in the s-c plane. An unsuppressed coupling c ~ 0.1 is allowed by the data, 
leading to a noticeable scale-dependence of the spectral index. The fit with 
c = 0.1 is less good than with a scale- independent spectral index, but still 
acceptable. 
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Figure 1: The plots show nominal 1- and 2-a bounds on n. In (a) the fit is 
performed fixing the reionization epoch zr, while in (b) is fixed instead the 
fraction / of matter which has collapsed at the epoch of reionization (the 
corresponding reionization redshifts, at best fit, are in the range 10 to 26) 
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Figure 2: The parameter space for the running mass model. As discussed 
in jH]] the model comes in four versions, corresponding to the four quadrants 
of the parameter space. In panel (a), we show the region allowed by observa- 
tion, in the case that reionization occurs when / ~ 1; the scale-dependence 
of the prediction for n is also displayed in this panel by the branches of the 
hyperbola 8sc = An = ns — ticobe, for the reference value An = 0.04. In 
panel (b), the straight lines correspond to ucobe = 1.2,1.0 and 0.8, and 
the shaded region is disfavoured on theoretical grounds 



